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We present our final results of the charmonium spectrum in quenched QCD on anisotropic lattices. Simulations 
are made with the plaquette gauge action and a tadpole improved clover quark action employing ^ = as I at = 3. 
We calculate the spectrum of S- and P-states and their excitation, and study the scaling behavior of mass splittings. 
Comparison is made with the experiment and previous lattice results. The issue of hyperfine splitting for different 
choices of the clover coefficients obtained by Klassen is discussed. 



1. Introduction 

Standard lattice QCD actions on space-time 
isotropic lattices encounter serious obstacles for 
heavy quarks with currently accessible lattice 
spacings because mass-dependent 0(ma) dis- 
cretization errors are very large. Aiming to re- 
duce such errors, Klassen[|l|,|| has proposed to 
employ anisoropic lattices with mat ^ 1 for 
heavy quark simulations. In this paper, we sum- 
marize our final results of the quenched charmo- 
nium spectrum using the anisotropic method . 
We also address the problem with hyperfine split- 
ting!^ that different choices of clover coefficients 
lead to disagreeing results in the continuum limit. 

2. Simulations 

We use the standard anisotropic gauge action 
given by Sg = ^Ei^/^oPss' + Co-P.*)- The bare 
anisotropy is tuned to obtain a desired value of 
the renormalized anisotropy ^ = as /at, adopting 
Klassen's parametrizationj|]. 
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Table 1 



Simulation parameters, a^ is fixed by ro = 0.5 fm. 
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For quark we use an anisotropic clover quark 
action: 

-Ks[M^ - ii)U,,xi^^^i + i'x+ii^ + liWlxi^x] 

+iKs[cs-il^xa-tjFij{x)ipx + Ct-il^xa-oiFoi{x)^px]}- (1) 

The bare quark mass is given by mg = \/2Kt ~ 
3/C — 1 with Q = Kt/Ks- For C, we adopt the tree 
level tadpole improved value for massive quarks. 
For clover coefficients Cs and Cf, we employ the 
values in the massless limit. We note that our 
choice of Cs is still correct for massive quarks 
because it has no mass dependence at the tree 
level [Q. The tadpole factors are determined as 
{Us) = {Pss')^^^ with Pss' the spatial plaquette 
and {Ut) = l. 

Simulation parameters are summarized in Ta- 
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Figure 1. cc spectrum with IP — IS input. 

ble 1. We adopt lattices with ^ = 3 and Lag ^ 
1.6 fm. Runs are made at four values of /? which 
correspond to Qs = 0.07-0.20 fm. For each /3, we 
measure S- and P-state meson correlation func- 
tions at two values of bare quark mass. Results 
are then inter(extra)polated to the charm quark 
mass where IS* mass has its experimental value. 
The lattice scale is set by either the Sommer scale 
To = 0.5 fm, IP- IS* splitting or 25-1^ splitting. 

3. Results 

In Fig.Q, we show results of the charmonium 
spectrum with the scale from the IP — IS' split- 
ting. Gross features of the spectrum are consis- 
tent with the experiment, e.g. splittings between 
Xc states are well resolved with correct ordering. 
The deviation of 2S masses from the experiment 
is in part ascribed to the quenching effect and in 
part to contaminations from higher excited states. 

3.1. Hyperfine splitting 

In Fig.^, we plot by filled symbols the lat- 
tice spacing dependence of the hyperfine splitting 
AM(l^S'i - l^So) for three inputs for the scale. 
Data at finite are extrapolated to the con- 
tinuum limit adopting an -linear ansatz. The 
results largely depend on scale inputs, and are 
much smaller than the experimental value {e.g., 
by about 30% with IP— IS input). Thus quench- 
ing effects are very large for the hyperfine split- 
ting. 

In the same figure, we also plot results by 
Klassen (open diamonds; f = 3)[^ and Chen 
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Figure 2. Hyperfine splitting. 
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Figure 3. Fine structure. 



(open triangles; ^ = 2)1^ with the same action. 
Their simulations differ from ours in that we de- 
termine the tadpole factor uq from the plaquette 
average and adopt for the parameter ( the tree- 
level tadpole improved value C"^^, while they use 
the mean link in the Landau gauge for uq and a 
non-perturbative estimate determined from 
the meson dispersion relation. Nonetheless, their 
results and ours, using the same scale rg, all con- 
verge to a consistent value of about 70 MeV in 
the continuum limit. 

3.2. Fine structure 

Figure ^ shows results of the fine structure 
AM{l^Pi - l^Po). The deviation from the ex- 
perimental value is smaller than that for the hy- 
perfine splitting (about 20% with IP — IS input). 
Our result with vq input is again consistent with 
those of Refs. fii. 



4. Effect of Cs for hyperfine splitting 

The results described so far all use the tad- 
pole improved value Cg = 1 for the spatial clover 
coefficient. In Refs.[Q,||, Klassen employed a dif- 
ferent choice Cg — ^Iv iv = ^o/O- He obtained 
HFS(a^ = 0,ro) « 90 MeV for the continuum 
limit of the hyperfine splitting, which is much 
larger than the result above HFS(as = 0, rg) « 
70 MeV with c~ = 1. We note that = l/i^ is 
correct only in the massless limit, while Cs = 1 is 
valid for any quark mass, at the tree level. 

To resolve this problem, we attempt an ef- 
fective analysis. The potential model predicts 
that the hyperfine splitting is due to the spin- 
spin interaction of quarks, which originates from 
the S • B term in the nonrelativistic Hamiltonian 
H^^. We therefore define a "tree- level effective 
hyperfine splitting" 
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is the tree level coefficient of the S • B term in 
H^^. The pole mass at Mi = log(l -I- mp) is 
inserted to normalize to unity in the continuum 
limit, and tildes denote the tadpole improvement. 

In Fig.| we compare the scaling behavior of 
HFS''^ (left panel) and the actual data HFS (right 
panel) for Cs = A similar comparison for 

Cs = 1 is made in Fig. ||. We find that results of 
HFS are qualitatively well reproduced by those 
of HFS'^'^. For c~ = l/v, HFS°* remains large 
even at (osMi)^ ~ 1, which suggests that the 
actual HFS should rapidly decrease as Cs ^ 0, 
and hence a naive estimation « 90 MeV[|l],|j from 
an a^-linear continuum fit is misleading for this 
case. On the other hand, HFS''^ is already close 
to unity for (asMi)^<l for Cs = 1. Thus an a^- 
linear continuum estimation (w 70 MeV) for this 
case appears much more reliable than that for 
Cs = l/v. 

5. Conclusions 

We have computed the charmonium spectrum 
accurately using quenched anisotropic lattices 
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Figure 4. HFS"'^ and HFS for c~ = l/v. 
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Figure 5. HFS°** and HFS for c~ = 1. 

with ^ = 3. We find that the spin splittings 
largely depend on the scale input and are smaller 
than the experimental values. Our results are 
consistent with previous results |^,|| when the 
same clover coefficients are used. We have also 
shown that a large hyperfine splitting reported 
in Ref. with a different choice of the clover 
coefficients is likely an overestimate arising from 
the continuum extrapolation. 
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